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TRANSCENDENTAL TRACE FORMULAS 
FOR FINITE-GAP POTENTIALS 

YURII V. BREZHNEV 



Abstract. We show that formulas differing from classical analogues of rational trace 
formulas for algebraic-geometric potentials occur in the theory of finite-gap integration 
of spectral equations. The new formulas contain transcendental modular functions and 
hypergeometric series. They result in transcendental relations for theta functions. 
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1. Intorduction 



Trace identities and formulas are known in the integration theory of spectral problems. 
This terminology originates from papers by Krein, Gel'fand, and Dickey on the spectral 
theory of the ordinary differential equation of the form tp xx — u(x)ip = Xip dating back 
to the 1950s. In 1975, stimulated by the appearance of finite-gap integration methods 
(in the spectral formulation) , V. Matveev established a new and important interpretation 
of these formulas. These relations, together with their t-isospectral modifications, are 
differential identities for exactly solvable potentials u{x). In particular, using one of the 
series of these formulas, namely, the famous formula 

u(x) = 2{ 7l (x) + ---+7 9 (x)} + const, (1) 

allows reconstructing the potentials from the quantities 7 fc (x), which satisfy integrable 
ordinary differential equations that are the Dubrovin equations. All the remaining trace 
identities follow from this "master" formula. Here, we interpret these relations or, to be 
precise, their analogues for higher-order spectral equations as formulas that express an 
exactly solvable potential in terms of the quantities J k (x). 
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It is known that a potential u{x) determines a Liouville-integrable finite-dimensional 
Hamiltonian system in which x plays the role of time. We pass to the canonical variables 
(p, q) by combining the potential and its derivatives: p = p(u, u x , . . .), q = q(u, u x , . . .) 
[8] . The subsequent transition to separation variables is the transition to the variables 7 fc 
and n k subject to the algebraic relation W(j k , fj, k ) = 0, which is called the spectral curve 
equation. Formulas of the inverse transition p = p{"f kl H k ), q = q( r f k , fi k ) and hence the 
formula of the form u = u(p, q) are well known [2] in Hamiltonian system theory, but we 
can clearly identify formula (1) with a differential analogue of an operator trace only in 
the Schrodingcr equation case. A "trace" interpretation of all other algebraic-geometric 
(finite-gap) operators is lacking. This is because both the spectral curve coordinates 7 and 
fi must participate in a general theory of variable separation because Abclian functions, 
generally speaking, are expressed in terms of both these variables. Second-order equations 
are the only particular case where u = 11(7). Only isolated examples are known for more 
complex spectral problems [4] , while a general scheme for constructing the principal trace 
formula 

u = R{ lk ,ii k ) (2) 
with a rational function R is still lacking, to the best of our knowledge. 

We devote this paper not to this problem but to the observation that in addition to 
representing potentials in the form of Abelian functions as in (2), we have other rep- 
resentations. They are provided by transcendental, not rational but still single-valued, 
functions of symmetric combinations of the parameters (7/-, The example we use to 
demonstrate this is interesting not only in itself but also because it leads to new tran- 
scendental identities for 0-functions. We recall that polynomial identities for O-functions 
are well known and are in fact uniformizing representations for Abclian varieties [3] and 
0-functional representations for Abelian functions (2) respectively. 

2. The rational trace formula 
We consider the third-order spectral problem 

*'" + u(i)*'-^'(a;)* = A«t. (3) 

It is known that if we associate this problem with a nonlinear integrable equation, then 
its isospectral deformations are described by the Kaup-Kupershmidt equation [11] 

25 

u t = u xxxxx + 5uu xxx + —u x u xx + 5u 2 u x . (4) 

In what follows, it suffices to consider only stationary solutions or, equivalently, the 
operator pencil 

A(A;[u])tf" + B(A;[u])tf' + C(A;[u])tf =/itf, (5) 
which depends linearly on A and commutes with (3). 

2.1. The algebraic curve. The compatibility condition of Eq. (5) with the spectral 
problem, i.e., with problem (3), customarily yields expressions for A, B, and C: 

A = -9A, B = u 2 + -u" + c, C = —A" + -Au- B' , 

2 3 3 

The constant is related to the stationary dynamics u — u(x — ct), and we have 
-9A*"+ (u 2 + ^u" + c)*' - (6u\ + 2uu" + t^'")* = M*- 
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The quantities A and fi are algebraically dependent, i.e., they belong to the algebraic 
curve 

[i 3 + Q(X;[u]) f i + T(X;[u]) = 0. (6) 

We obtain the equation of this curve by eliminating ^ from (3) and (5), with the result 

/ 9 \ 27 1 

fi 3 + [27 c\ 2 --KX + E 2 /i + 729A 5 + £iA 3 + — uK X 2 + E 3 X+ -(u" + 2u 2 + 2c) 2 K = 0. 
V 2 / 2 8 

We present expressions for E^ below and here note that not only integrals of stationary 

equation (4) but also the equation itself enter the spectral curve equation: 

OK 

K = cu' + u (V) + 5uu"' + y uV + hu2u '- 

By virtue of (4) the constant K must vanish identically. But the expression u" +2u 2 +2c = 
const also does not contradict it. If the constant in the right-hand side of this equation 
is nonzero, then the compatibility with the equation K = results in the trivial answer 
u = const. Otherwise, we have the degenerate curve 

^ + 27cX 2 n + 729 A 5 - 27 (^u' 2 + u 3 + Zcuj A 3 = 0. 

This curve has the genus 3 = 1, the corresponding solution u = — 3p(x + f^) can be 
easily calculated, and we are not interested in it. In the nondegenerate case, curve (6) 
becomes 

fi 3 + (27 cX 2 +E 2 )fi + 729 A 5 + 81 £i A 3 + E 3 X = 0, (7) 

it has the genus g = 4, and the problem of determining the corresponding potential is 
nontrivial. 

2.2. The rational trace formula. We first obtain the standard trace formula, i. e., the 
representation of the potential in form (2). Let us determine the quantities j k . In the 
finite-gap potential theory, these objects are the poles of ^(x), which arc independent of 
A. Eliminating the ^-function and its first derivative from Eqs. (3) and (5), we obtain 
the expression 

Here Q(X) := 27c A 2 + E\ ([«]), and the differential polynomial II = II(A; [u]) is 

3- 7 II(A; [«]) = A 4 - -u'X 3 + —(u (IV) + 5uu" + Au' 2 + 2u 3 + 2cu)A 2 + 
6 18 

+ ^rA u '" + 4uu')(u" + 2u 2 + 2c) A + ^-{u" + 2u 2 + 2c) 3 . 
324 v M ' 3 6 8 

We now present the Novikov integrals Ek'. 

E x = u {lv) + 5uu" + ^ u' 2 + ~ u 3 + cu, (10) 

E 2 = —u"'(u"' + 8uu') - -uu" 2 + —u" (Ex + 'iu' 2 - 2{)u 3 - 12 cu) + ■■■ , 
4 4 24 

E 3 = --u"'u'(u" + 2u 2 + 2c) + -u" 3 + - (u 2 + c)u" 2 - 3u' 2 u" + ■ ■ ■ , 
4 4 4 

where we truncate the expressions for E 2 , E3 because we do not need them in what 
follows. 



0) 
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The separation variables j k determine the poles of logarithmic derivative (8), i. e., they 
factor the polynomial II. We therefore define 

n = 3 7 (A - 7l )(A - 72 )(A - 73 )(A - 7 4 )- (11) 

If we define the second coordinate [i k = fJ- k (x) by the formula 

Vl + Q{l k (x))v k +T{ lk (x)) = 0, 

then we obtain an analogue of the Dubrovin equations from (8) in the case under consid- 
eration. Indeed, passing to the limit A — > j k in (8), we obtain 

dlk _ Ik 3/4 + Qilu) 



dx 81 n(7 fc "7,) 



(12) 



Additional explanations can be found in [4] , where the general method can also be found 
for obtaining such equations and the second coordinate of the zero of the polynomial II: 

K + 2u 2 + 2c) 2 

M ^ 3M7fe + 36^ ■ (13) 

Comparing (9) and (11), for instance, we obtain the identities 

4 

w' = 6^7 fe , u (lv) +5W + 4u' 2 + 2u 3 + 2cu = 18^ 7fc7j . (14) 

l k>j 

Composing symmetric combinations of Eqs. (12) and using formula (13) for the second 
coordinate, we can eliminate all the derivatives of the potential and thus obtain a linear 
equation for u. This is the desired trace formula. Omitting the details, we eventually 
obtain 

1 &k 

The last simplification of Eqs. (12) is to reduce them to quadratures. Again composing 
symmetric combinations of the right-hand sides of (12), we obtain the equations 

V ^ = V lkdlk =o 

^3/4 + Q( 7fc ) ' V fc3 ^ + W ' 

E H k dlk n >p ijdlk - — rl 

1 *3A*1 + Q( 7fc ) " ' ^3ti + Q( lk )~-81 X > 

which, as expected, result in a Jacobi problem, i.e., in the problem of inverting four 
holomorphic integrals on curve (7): 

dx A r \d\ 



k / 3,,2 + om =ai ' 



/* J 3m 2 + Q(A) ^ 7 3^ 2 + Q(A) 



'2- 



a (7*. 4 (7fciMh) 

A h /idA A / A 2 dA 1 

^ 7 3 M 2 + Q(A) ~ fl3 ' ^ y vTow " a4_ 8i x - 

In what follows, we consider the dependence 7 fc = 7fc (a;) to be defined from these equa- 
tions. 
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Wc note that this scheme for obtaining the Dubrovin equation and trace formula from 
the definition of the differential polynomial II and "J is intrinsic for general finite-gap 
operators. The definitions of the quantities j k , the Dubrovin equations, and various 
relations for the potential and its derivatives (different versions of the trace formula) are 
"simultaneously entangled" in the general set of determining polynomial relations. 

3. The transcendental trace formula 

Deriving the preceding formulas, we manipulated with polynomials, which is common in 
polynomial ideal theory [7] . The bases of such ideals are not unique, which means that we 
can presumably find other relations determining the potential. This is indeed the case. 
Eliminating the higher derivative u^ IV) from the Novikov integral E\ and from the second 
identity in (10), we obtain the quadratic equation for u': 

3u' 2 + Au 3 + 12cu + Ex = 216 ^ . 

k>j 

But v! can be expressed uniquely in terms of 7 by virtue of the first identity in (14). We 
then obtain the cubic polynomial in the variable u: 

4 

u 3 + 3cu + 3E 1 + 27Y^ ll = 0. (16) 
1 

It can be interpreted as an algebraic (implicit) variant of the principal trace formula. 

On the other hand, finite-gap potentials are single- valued functions of the variable x. 
They satisfy autonomous ordinary differential equations, pass the Painleve test, and are 
expressed in terms of theta functions. It is also known (but seldom used) that roots of any 
3rd or 4th degree polynomial can be expressed analytically via its coefficients in terms 
of elliptic functions [1]. For exhaustive information on this topic see [5]. In our case the 
situation simplifies since we may directly consider polynomial (16) as a Weierstrass cubic 
polynomial 

4 

Au 3 + Ylcu + Y1E X + 108 V jlix) = 4m 3 - au - b = 

— ' k 
1 

= 4(u-e)(u-e')(u-e"), 

where the points e, e', and e" depend on the coefficients a and b, which in turn depend 
on x. The roots of this polynomial are 

u = {e, e', e"} = p(u k (x); a, b) , 

where 

a = -12c, & = -108X)7fc(y)-L2£i, 
are expressed in terms of the theta constants [1] 

e= ^K(r)+*i(T)}, 

e' = -i^K(r)+^(r)}. (17) 
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and the quantities uj and t = ui' /ui arc expressed in terms of the coefficients a and b. ft 
suffices to take just one such expression and express the ^-constants in (f7) using the 
classic series [f]: 

oo oo oo 

Mr) ■= e^ T J2 k e ik2+k)niT , Mr) := E^™' Mr) ■= ^(-l)*^. 

— oo — oo — oo 

(we let the symbol c z denote the exponential function to distinguish it from the Weier- 
strass e-points). 

Wc thus obtain a classic elliptic modular function inversion problem, i. e., the problem 
of expressing the periods (2ui, 2ui') of the elliptic curve w 2 = 4z 3 — az — b in terms of its 
coefficients. The solution scheme is known. We must find the root r of the transcendental 
equation 

J(r) = - 



3 - 27b 2 ' 

where J(r) is the classic modular Klein function [f, 10]. After finding the root r, we 
calculate the half-periods (uj, uj') by the formula 



2 a 9zi T ) 

uj = T—j—r, uj = ra, (18) 
o .92 ( T ) 

where g-zsir) are the known modular forms [10]. We have the Eisenstein series, the 
Hurwitz-Lambcrt series, or the ^-constant representations for these forms: 



S 2 (T) = ^{tf!(T) + tf!(T)+tf!(T)} 

ffa(r) = ^ {^(r) + ^(r)}{^(r) + ^(r)}{^(r) - Mr)} 



(19) 



It remains to write the expression for the root r. Strangely enough, an explicit analytic 
solution of this classic problem is lacking 1 , although it is common knowledge that this 
solution can be written in terms of ratios of hypergeometric 2-F i-series. Because the series 
2-F1 (J) converges only inside the unit circle, the formulas differ depending on whether 
\J\ > 1 or \J\ < 1. Altogether, this results in cumbersome expressions also containing 
functional series in the logarithmic derivative of the Euler T-function. For example, the 
corresponding expression (22)-(27) in Sec. 14.6.2 of [10] is half a page long (even if we 
forgive the incorrectness of formula (23)). But the problem admits a simple solution. 

We use the well-known fact that the function J is closely related to the hypergeometric 
equation of the form 

J(J-lW + l (7J- 4)^ + ^ = 0. 

Solutions of general hypergeometric equations are hypergeometric functions 2^1(0, b; c\z), 
but these functions can be rewritten in terms of classic special functions under special 
restrictions imposed on the parameters a, 6, and c. We mean not reductions of the 
function 2-F1 but the cases where the 2-Fi-series admits a quadratic transformation, and the 
hypergeometric equation then reduces to a two-parameter equation, e. g., to the Legendre 
equation. This is true for the above equation. Its solution is a linear combination 

x/j = Vj{apz(VT~j) + bq^v 7 ! - 7)} 



^Unfortunately, solutions of this problem are given incorrectly both on page 789 of encyclopedia [6] 
and in the end of Sect. 11 of the remarkable book [!]; this is not due to misprints. 
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of the Legendre functions with the parameters (y, fi) = (—5, 3) (see [9] for detailed 
information about these functions). From the above, we conclude that the formula 



aP(VT~J) + bQ(VT-J) 

cP(vT^7) + dQ(VT^~J) 

holds at certain numerical values of {a, 6, c, e£} (we omit the indices is, \i of the Legendre 
functions for brevity). Wc present only the final answer; see [5] for derivation. 

Proposition 1. For the elliptic curve w 2 = 4z 3 — az — b, the period ratio r = cj'/cj is 



Q(v / T^J) j ' a 3 -276 2 

Remark 1. If J is substituted for J(r) in equality (21), then this equality becomes an 
identity that holds for all r in the upper half-plane H + . The left- and right-hand sides of 
this identity must be calculated up to the action of the total modular group PSL^Z) =: 
r(l). Relation (21) is therefore a r(l)-equivalent of the solution of the modular inversion 
problem in the Legendre representation. Namely, given the elliptic curve in the Legendre 
form w 2 = (1 — z 2 )(l — k 2 z 2 ), we can calculate the elliptic modulus r by the celebrated 
Jacobi formula 

t=[ iM m ° dr(2) ' (22) 

where K and K' are the complete Legendre elliptic integrals of the first kind [1, 10, 13]. 

Returning to the potential, we use formulas (17) and (18) with (19) taken into account 
to obtain the desired transcendental "trace formula" : 

«(*) - ^ ^ + ^ 7fc (x)| m + n{T)m(T) _ n{T)] , (23) 

where r is the function 



V Q(vT^J) y 4 C 3 + 9{i? 1 +9E7 fe 2 (x)} 2 

and we let the symbol T denote the operation of setting a number to the fundamental 
domain of the group. Formulas (23)-(24) do not contain the coordinates /i fe . Verifying 
that function (23) indeed satisfies relation (10) or a stationary version of Eq. (4) is a good 
and rather nontrivial exercise. 

4. Transcendental identities for 6-functions 

The results in the preceding section demonstrate that we can analogously obtain many 
more examples of Abelian functions that are radicals or roots of algebraic equations. For 
example, the last summand in polynomial (9) provides a single-valued representation of 
the Abelian function radical of the form 

18^7i7 2 7 3 74 =u" + 2u 2 + 2c, 

and the general mechanism for constructing "transcendental traces" is based on eliminat- 
ing derivatives of the potential from the total set of relations determining the integrals 
and the variables j k . Solving the obtained equations in single- valued functions (which are 



8 



YU. BREZHNEV 



always theta functions), we obtain the answer. These equations (and even possibly their 
orders) are not uniquely defined and may depend on the genus, but transitions between 
them are just different representations of the same Abelian function. They differ even in 
the (7, /^-representation because the coordinates 7 and /z are algebraically related. 

Now equating the transcendental and rational representations, we obtain an unusual 
identity containing theta functions. Indeed, comparing Eqs. (15) and (23), for example, 
(23) 



3 Z^; 



llH 3 f * \ tff(T)+tf§(T)+tfj(T) 



n(7 fc -7,) 2c\ v i IWw+^DWw-^t)}' 



we can in principle represent the Abelian functions rational in (7, fx) in this formula in 
the form of 8-function relations [3] . The obtained relation for 6-functions of genus g = 4 
is a transcendental identity containing the theta constants $(t(x)), although curve (7), 
generally speaking, is not related to any elliptic curve. In addition to the curve moduli, 
this identity also contains an arbitrary parameter x. We can also move all symmetric 
functions that are rational in (7, /1) to the left-hand side of the equality. We then obtain 
an Abelian function for the Jacobian of trigonal curve (7) in the left-hand side and an 
"exotic" representation of this curve in terms of the elliptic theta constants of the Legendre 
functions in the right-hand side: 

u 2 = -2c x 2V ' , 3V 8 J .\ Is/ ■ ~ 3c - 26 

In turn, arguments of the Legendre functions (hypergeometric functions) in this repre- 
sentation are expressions containing another Abelian function via the quantity yl — J, 
i.e., J2^k( x )> by virtue of formula (24). This function is represented by its own rational 
fraction of 0-functions of the type Q(xU + D). The nature of transcendental relations 
is therefore rather involved and still needs to be understood. The transcendency is pre- 
served even when the Jacobian of curve (7) splits completely into elliptic curves, and the 
multidimensional 0-functions become the elliptic ^-functions. Polynomials of type (16) 
then remain nondegenerate, and it is natural to expect the appearance of transcendental 
^-identities. 

The type of transcendency obtained above is related to the elliptic modular inversion 
because this inversion is generated by a solution of an algebraic equation with application 
of the elliptic modular functions. The general theta functions and constants appear in 
more-general methods (see, e.g., Umemura's contribution to [12]), and it is therefore 
natural to expect the appearance of other types of transcendental identities. 

In conclusion, we note that the simultaneous and transcendental unification of theta 
functions and hypergeometric scries is in fact not unexpected, because a natural feature 
of 2-F1 series are their monodromy groups, which are in turn automorphisms of analytic 
automorphic functions. All such functions that are currently known are ratios of theta 
constants. Informally speaking, the "theta" geometry and "hypergeometry" are mutu- 
ally inverse (they realize inversions of fractions of type (20)), and our transcendental 
identities reflect both this fact and the properties of functions of the type Q(xU + D). 
Namely, Abelian manifolds are parameterized as algebraic manifolds by polynomials in 
theta functions while linear sections of Jacobians appear from integrable equations. A 
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partial illustration is already provided by Jacobi equation (22) rewritten in the form 



Considering this equation in a sufficiently small vicinity of a point r , we can drop out 
a "not completely analytic" operation of bringing a point into the fundamental domain 
of the group, and this equation then becomes an exact analytic transcendental identity. 
We can proceed in the same way with the theta versions of our relations (25) and (26) 
interpreting them as functions of x. 
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